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Abstract: Due to the characteristics of high precision for a long duration, the rotary 
strapdown inertial navigation system (RSINS) has been widely used in submarines and 
surface ships. Nowadays, the core technology, the rotating scheme, has been studied by 
numerous researchers. It is well known that as one of the key technologies, the rotating 
angular rate seriously influences the effectiveness of the error modulating. In order to 
design the optimal rotating angular rate of the RSINS, the relationship between the rotating 
angular rate and the velocity error of the RSINS was analyzed in detail based on the Laplace 
transform and the inverse Laplace transform in this paper. The analysis results showed that 
the velocity error of the RSINS depends on not only the sensor error, but also the rotating 
angular rate. In order to minimize the velocity error, the rotating angular rate of the RSINS 
should match the sensor error. One optimal design method for the rotating rate of the RSINS 
was also proposed in this paper. Simulation and experimental results verified the validity and 
superiority of this optimal design method for the rotating rate of the RSINS. 

Keywords: rotary strapdown inertial navigation system; rotation modulation; rotating 
angular rate; optimal design 
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1. Introduction 

The operation of aninertial navigation system (INS) depends on Newton's second law [1]. The inertial 
navigation is the process whereby the measurements of inertial sensors (gyroscopes and accelerometers) 
calculate the position and velocity of the vehicle. Unlike other types of navigation systems, the INS is 
an entirely independent navigation system, since it is not dependent on the transmission of signals from 
the vehicle or reception from an external source. However, the fatal flaw of the INS is that the navigation 
error caused by the bias of the inertial sensors increases with time. Therefore, the precision of the INS 
is subject to the precision of inertial sensors [2,3]. 

For the original INS, the inertial sensors are mounted on a stable platform, which is mechanically 
isolated from the rotational motion of the vehicle [4] . Due to the mechanical complexity of the platform 
system, a new type of INS, which is called the strapdown inertial navigation system (SINS), is replacing 
the traditional platform inertial navigation system. The SINS has the advantages of lower cost, smaller 
size and greater reliability compared with platform systems, since the inertial sensors are attached rigidly 
to the body of the vehicle [5]. 

Since the navigation error caused by the bias of the inertial sensors increases with time, researchers 
proposed the rotation modulation method to reduce the impact on the SINS, which is installed on the 
vehicle through the turntable. The rotation modulation method modulates the bias of inertial sensors into 
a zero-mean and periodical form to inhibit the error propagation, and then, the localization precision of 
the SINS will be enhanced effectively. This kind of rotary strapdown inertial navigation system (RSINS) 
has been applied by the U.S. Navy since the 1960s [6-8]. However, the disadvantage of the rotation 
modulation is that it cannot modulate the installation error and the symmetric scale factor error of the 
inertial sensors, but it also produces instantaneous velocity error [8-10]. As a result, the design of the 
rotating angular rate has became a hotspot of rotation modulation. 

The RSINS accuracy influenced by the non-uniform rotation angular rate is analyzed by Dr. Che [11]. 
The rotating scheme design method of the non-uniform rotating angular rate is proposed in [12]. 
However, until now, no one has been able to figure out an effective optimal rotating angular rate design 
method of the RSINS in theory. Therefore, in this article, a novel rotating angular rate design discipline, 
which will fill the gap, was proposed. We focused on the optimal design of the rotating angular rate 
utilizing the Laplace transform and the inverse Laplace transform in this article, and this paper is 
organized as follows. In Section 2, the principle of the rotation modulation and the inertial sensor 
error model are introduced. The design method of the rotating angular rate is proposed in Section 3. 
Simulations and experiments are carried out in Section 4 and 5, respectively, to verify the correctness of 
the design method. 

2. Principle of the Rotation Modulation and the Error Model 

2.1. Rotary Strapdown Inertial Navigation System 

As the SINS is installed on the vehicle through the turntable, the structure of the RSINS, which is 
shown as Figure 1, is different from the SINS. Since the middle of the last century, rotation modulation 
has been widely applied to the SINS [9,10]. 



Sensors 2014, 14 



7158 



Classified by the number of the rotating axes, the RSINS can be divided into single-axis RSINS and 
the dual- axis RSINS [13]. However, the purpose is the same, i.e. to compensate for the effect of inertial 
sensor errors without external information [14]. 

Additionally, we must say that any technique has its limitations, and rotation modulation is no 
exception. One fatal shortcoming of rotation modulation is that it will produce instantaneous velocity 
error that cannot be ignored during the rotating period, and the error exists in any rotating mode [8,15]. 

Figure 1. Schematic of the rotary strapdown inertial navigation system (RSINS). 

I 

Navigation Calculation ■ Information Display 




Turn table 

7777777 



The essence of rotation modulation is the error self-compensation [16-18]. It modulates the inertial 
sensor errors in terms of periodical rotation. It is an effective method to enhance the navigation accuracy 
under low-accuracy sensors conditions [19-21]. The inertial sensor errors neither disappear nor change, 
but offset during the whole rotation period [21-24]. 

As the inertial measurement unit (IMU) connects with the vehicle indirectly through the turntable, 
a new frame, the sensor frame (s frame), is introduced in the RSINS. It is supposed that the body 
frame {h frame) and the navigation frame {n frame) are coherent. This means = I and = C^. 
The output of the inertial sensors is transformed from the s frame to the n frame by a direction cosine 
matrix. Suppose the s frame and n frame overlap at the initial time, and the IMU is controlled to rotate 
uniformly at a certain speed, cOr, along the vertical direction. The direction cosine matrix from the s 
frame to the n frame can be expressed by: 



nJ 



{CD 



s\T 



COS Urt 
sin UJr-t 

0 



— sin ojrt 0 

cos UJrt 0 

0 1 



(1) 



The modulated constant bias of the gyroscope in the n frame [cei ^ni ^uiY tie described by 
the following equation. 













= c: 















e^O cos Urt — Eyo sin COrt 
ExO sin Urt + Eyo COS OJrt 



(2) 



where e^o, £yo^ ^zo denote the constant bias of the three orthogonal gyroscopes. From Equation (2), 
it is obviously known that the horizontal constant bias of the gyroscope can be modulated into periodic 
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and zero-mean signals by periodically rotating the IMU. Analogously, it has the same result for the 
accelerometer. This technique can be utilized to modulate the inertial sensor error, which is vertical to 
the rotating axis, in order to enhance the positioning accuracy of the RSINS. 







cos Ulrt 


— sin uirt 








Va;0 cos Ulrt — VyO Sin Ulj-t 


Viv 




sin uirt 


cos Ulrt 








VyQ sin uij-t + V xQ cos uij-t 



2.2. The Sensor Error Model of the RSINS 

The inertial sensor error mainly includes the constant bias, the scale factor error and the installation 
error of the optical gyroscope and the constant bias of the accelerometer. In Section 2, we will focus on 
the relationship between the inertial sensor error and the performance of the RSINS. 

2.2.1. Constant Bias of Optical Gyroscopes and Accelerometers 

The modulated constant bias of the gyroscope and the accelerometer in the n frame can be described 
by Equations (2) and (3), respectively. 

(3) 

L Lklrl -\- V xO cos ljjj.t 

where Va;o and V^o are the constant bias of horizontal accelerometers. 

2.2.2. Scale Factor Error of Optical Gyroscopes 

The equivalent bias of gyroscopes, e', caused by the scale factor asymmetric error of gyroscopes can 
be described by: 



(4) 



where Kgx,Kgy and Kg^ denote the scale factor asymmetry errors of gyroscopes along the x, y and z 
axes, respectively; denotes the Earth's rotation rate; L denotes the local latitude; [Ux cOy denotes 
the output of gyroscopes, and its theoretical value is: 
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(5) 



Therefore, the modulated equivalent bias of gyroscopes, e', in the n frame can be described by: 
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(6) 



2.2.3. Installation Error of Optical Gyroscopes 

The inertial sensors should be installed orthogonally in theory. However, the actually s frame is not 
orthogonal and cannot coincide with the h frame, as shown in Figure 2. The small-angle 5ij is defined 
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as the installation error between the i axis of the s frame and the j axis of the h frame. Therefore, the 
installation error matrix can be described as: 
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(8) 



Therefore, the modulated equivalent bias of gyroscopes, e", in the n frame can be described by: 
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(9) 
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2.2.4. Error Model of Inertial Sensors in the RSINS 

Concerning the study above, the error model of the RSINS can be described as 









Vat 




Vat 


0 
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£e 




^El +^E2 +£E3 
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£U1 + ^U2 + £(73 



Wi{t) + W2{t) + W3{t) 



(10) 



where Wi{t), W2{t) and W3{t) represent the constant bias error of inertial sensors, the scale factor error 
and the installation error of gyroscopes, respectively. 



Vxo cos ojrt — Vyo sin OJrt 

VsO sin UJrt + VyO COS OUrt 

0 

£xO COS OJrt — £yO Sin Uj-t 
£yO COS OJrt + £xO Sin OJrt 
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0 
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W^it) 



0 
0 
0 

e'I cos OJrt — e'y sin ijj.,t 

£'' cos OJrt + £'L sin Urt 



It can be seen from Equation (10) that the error model of the RSINS is a function that depends on not 
only the error of inertial sensors, but also the rotating angular rate. Therefore, the performance of the 
RSINS will be influenced by the rotating angular rate. 



3. Optimal Design Method for the Rotating Angular Rate in the RSINS 

The previous analysis shows that the error model of the RSINS bears on both the inertial sensor error 
and the rotating angular rate. Therefore, this section will investigate the relationship between the rotating 
angular rate and the RSINS performance. 

3.1. The Velocity Error Equation of the RSINS 

Based on the spherical Earth model, ignoring the movement of the vehicle, suppose the IMU rotates 
along the vertical axis at ojr', the error model of the RSINS can be described as: 



X{t) = AX{t) + Wi{t) + W2{t) + Ws{t) 



(11) 
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where X(t)=[(5VB 5Vn 5L i^e v^c/]^ and: 
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where 5Ve and 5Vn represent the velocity errors in the east and north, respectively; SL represents the 
error of the local latitude; ifE, and ipu represent the misalignment angle of the east, north and up; R 
represents the radius of the Earth; g represents the gravitational acceleration. After the Laplace transform 
and ignoring the initial value. Equation (11) is transformed to the following format: 



X{s) = {si - A)-' [Wi{s) + W2{s) + Ws{s)] 



(12) 
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As the values of rrij^^j = 1, 2, ■ • ■ , 6) do not affect X{s), the elements of the third column are 
represented with '*'. Wi{s), the expression of Wi{t) in the frequency domain, can be expressed as: 



W^{s) 



where: 



LOr'VxO + s'VyO 
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Therefore, the velocity error of the east and north in the frequency domain can be expressed as: 



(13) 



SVe{s) = Ml • Wi{s) + Ml ■ W2{s) + Ml • ^^3(5) + Mi • ^4(5) 
SVNis) = M2 • Wi{s) + M2 • W2{s) + M2 • W^is) + M2 • W^s) 

where M^ • Wi{s){i = 1,2), M^ • W2{s){i = 1,2), M^ • W^3(s)(i = 1,2) and M^ • W4{s)(i = 1,2) 
represent the velocity error caused by the constant bias of gyroscopes, the velocity error caused by the 
constant bias of accelerometers, the velocity error caused by the gyroscope scale factor error and the 
velocity error caused by the gyroscope installation error, respectively. 



3.1.1. Velocity Error Caused by the Constant Bias of Gyroscopes 

The east velocity error caused by the gyroscope constant bias can be expressed in the frequency 
domain as: 

dV^is) = Ml • Wi{s) 

guf^ sin L ■ cos Le 



zO_ SgUie sinL(g£^o - Ur^yo) g(^ie Sm Ls{seyO + UJr^xo) '^^'^^ 

- S{S^ + ^2) (52 + + (,2 + ^2)(,2 + ^2J(,2 + ^2) (,2 + ^2)(,2 + ^2J(,2 + ^2) 

After the inverse Laplace transform and the simplification, we can get the east velocity error caused 
by the gyroscope constant bias in the time domain shown as: 



SV^ (t) =kii cos Ugt + A;i2 sin cust + kis cos cort + ku sin cori + ki5 cos coiet + kie 



(15) 
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AEliuJr) = \lkl^ + fc?2 + + + + (16) 

where: 

_ ^(jf^cosLsinL . _ 9 i, _ 9 i, _ 9 

Ku — 7 £zO', 1^12 — £zO', — :^£xO', Ku — ^^yO] 

gsinLcosL gsinL cos L 

kl5 = 5 £20; M6 = 5 £zO- 

The north velocity error caused by the gyroscope constant bias can be expressed in the frequency 
domain as: 

""^'^^ is' + U;'r){s'+ corns' +U:l) S{S' + U^S' + ^l) is' + CO^r){s' + CoDis' + Uj) ^ ^ 

After the inverse Laplace transform and simplification, one can get the north velocity error caused by 
the gyroscope constant bias in the time domain as: 

SV^{t) =k2l cos CUst + ^22 SmCUgt + /C23 COS CUrt + ^24 sin tU^t + k25 COS CUiet + ^26 sincjjg^ (18) 

The amplitude of the north velocity error is expressed as: 



Am M = \ + kl^ + kl^ + kl^ + A;2 + kl. (19) 



where: 



_ gujie sin L _ guie cos L _ 1 _ 1 

1^21 — 5 £xO'i 1^22 — ^ £zO; 1^22, — I^^xQ] i^2A — -^£yO] 

_ UiegsvuL _ gcosL 

K25 — 5 £xO] 1^26 — 5 £zO 

It is known that the amplitude of the velocity error caused by the gyroscope constant bias decreases 
with the increase of the rotating angular rate from Equations (16) and (19). Figure 3 shows the 
simulation results. 

3.1.2. Velocity Error Caused by the Constant Bias of Accelerometers 

The east velocity error caused by the accelerometer constant bias can be expressed in the frequency 
domain as: 

5Vi{s) = M, . W.is) = ^(;y.o - c..y,o) (20) 

After the inverse Laplace transform and simpUfication, we get the east velocity error caused by the 
accelerometer constant bias in the time domain and shown as: 



5 V| (t) = /cai cos CUst + k^2 sin <^st + k^^ cos Uj-t + ^34 sin 



(21) 
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Figure 3. The relationship curve of the amplitude of the velocity error, the constant bias of 
the gyroscope and the rotating angular rate. 




The amplitude of the east velocity error is expressed as: 

Ae2M = \l + kl^ + fcfs + (22) 

where: 

^31 = Vj,o; ^32 = ^Va;o; ^33 = V^o! ^34 = ^ yO 

iOr iOr 

The north velocity error caused by the accelerometer constant bias can be expressed in the frequency 
domain as: 

After the inverse Laplace transform and simplifying Equation (23), we get the east velocity error 
caused by the accelerometer constant bias in the time domain as: 

SV^ (t) =kii sin ojrt + ^42 cos Ust — ^43 cos Urt — ku sin cOst (24) 

The amplitude of the east velocity error is expressed as: 



^iV2(wr) = ykl^ + kl2 + kl^ + kl^ (25) 

where: 

/C41 — ; /v42 — ; /C43 — ; ^4 — 5 — 

UJr OOr OUr ijO^ 
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The conclusion is obtained from Equations (22) and (25) that the amplitude of the velocity error 
caused by the accelerometer constant bias decreases with the increase of the rotating angular rate. 
Figure 4 shows the simulation results. 

Figure 4. The relationship curve of the amplitude of the velocity error, the constant bias of 
the accelerometer and the rotating angular rate. 




3.1.3. Velocity Error Caused by the Scale Factor Error 

The east velocity error caused by the scale factor error can be expressed in the frequency domain as: 

5V|(s) = Ml ■ 1^3(5) 

gUJrOJie COS Li^S^UJie sin L — S^Ur — OJrOjff, cos^ L) 

^ {s^ + uI){s^ + ujI){s^ + uIY (26) 

s^gujie COS L(s^ + uJrOJie sin L + (jjf^ cos^ L) guf^ sin L cos L{uie sin L + Ur) 

After the inverse Laplace transform and simplifying Equation (26), we get the east velocity error 
caused by the scale factor error in the time domain as: 

5V|(t) =A;5i sin cort + k^2 sin Ust + /css sin Uiet + k^^t (27) 
The amplitude of the east velocity error is expressed as: 



AezM = ykli + kh + kh + Ha (28) 

where: 

guf^ sin L cos L guie cos L gcur sin L cos L 

^51 = E Kgz] K52 = n Kgx] ^53 = 5 ^Qz] 

gur sin L cos L 
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The north velocity error caused by the scale factor error can be expressed in the frequency domain as: 

5VIj{s) = M2-W^{s) 

guie COS LurS^s"^ + UrUJie sin L) s^gujie COS L{(jJr + ojie sin L) 

(s2 + a;2)(^2+^2J(^2+^2)2 9-' (^2 + ^2) (^2 + ^2J ^^2 + ^2)2 dV (29) 

(/Wje cos I/(a;ie sin L + ujr) 

After the inverse Laplace transform and the simplification, we get the east velocity error caused by 
the scale factor error in the time domain as: 

SV'^it) =kei cos COgt + kQ2 cos COrt + ^63 (30) 

The amplitude of the east velocity error is expressed as: 

AN3iuJr) = yjkl + kl^ + kls (31) 

where: 

_ gujfe'^r cos L _ gur COS L _ gur cos L 

fCei — 2 J^gz', /C62 — 5 J^gz', ^63 — 5 J^gz 

ujf 

From Equations (28) and (31), we can draw a conclusion that the amplitude of the velocity error 
caused by the accelerometer constant bias increases with the increase of the rotating angular rate. 
Additionally, we plot the relationship between the amphtude of the velocity error, the scale factor error 
and the rotating angular rate as Figure 5. 

Figure 5. The relationship curve of the amplitude of the velocity error, the scale factor error 
and the rotating angular rate. 
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3. 1.4. Velocity Error Caused by the Installation Error 

The east velocity error caused by the installation error can be expressed in the frequency domain as: 

[gujie sin Ls^ — gi^s"^ + uf^ cos^ L)ujr\ (cuie sin L + ujr) 
uiie COS Ls[guie sln Ls^ — g{s^ + uf^ cos^ L)u}r\ ^ 

ujie COS Lur [goOie sin Lscj^ + ^^(s^ + oofg cos^ L)s] (32) 

{ujie sin L + Ur) [gojie sin LscUr + + cufg cos^ L)s] 
+ .(52 + a;2)(52 + a;|J(.2 + ^2) 

5ia;fg sin L cos^ L g'O'fg sin L cos^ Loj^ 

After the inverse Laplace transform and simplification, we get the east velocity error caused by the 
installation error in the time domain as: 

SV^{t) —k-ji sin Urt + k'j2 cos Urt + /Cys cos Ugt + k-ji sin Uiet + A;75 cos Uiet + /Cye (33) 

The amplitude of the east velocity error is expressed as: 

Ae^M = \J k^i + k^2 + ^Is + ^74 + ^75 + ^76 (34) 



where: 



S"^j/a; , _ 5'<^a:2/ , _ fifWje sin L COS^ Lt^^j, 

/i^72 — ^; />^73 — 2 ; 



K'll — 5-, «^72 — 5-, «^73 

g sin LSxy 



_ g sin LSyx 

9 J '^75 — 9 ) '^'76 — q c 



A;, 



The north velocity error caused by the installation error can be expressed in the frequency domain as: 

g{s^ — UrOJie sinL)(a;je sin L + cUr) r. s^gojie cosL(s^ — UrOJie sin L) 

{ur + oJie sin L)s'^guJie COS Lur r. gs{ur + ooie sin LY 

" {s^ + UJ^,){s^ + coDis^ + a;2)2 + (s^ + a;2)(,2 + ^2 )(,2 + ^2) 

sguie^ cos^ L g^r^fe L 



(35) 



After the inverse Laplace transform and simplification, we get the north velocity error caused by the 
installation error in the time domain as: 

5Vj^{t) =ksi coscurt + k82 coscust + kss COS LOiet + k^A sincu^t • t (36) 
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The amplitude of the north velocity error is expressed as: 

ANi{uJr) = \l kl^ + fcia + Zeis + ^84 (37) 

where: 

Kgl — '''82 — 5-^, ^83 — n o , ^84 — 

From Equations (34) and (37), we can draw a conclusion that the amplitude of the velocity error 
caused by the installation error increase with the increase of rotating angular rate. Additionally, we plot 
the relationship between the amplitude of the velocity error, the installation error and the rotating angular 
rate as Figure 6. 

Figure 6. The relationship curve of the amplitude of the velocity error, the installation error 
and the rotating angular rate. 




3.2. Optimal Design Method for the Angular Rate 

Section 3.1 is focused on the relationship between the velocity error, three different kinds of inertial 
sensors errors and the rotating angular rate. In this part, we will propose the optimal design method of 
the rotating angular rate. 

When the RSINS rotates at a constant speed, the amplitude of the velocity error can be expressed as: 

A{Ur) = \lAl{uJr) + Al{Ur) + ^K^r) + ^l(c^r) (38) 

where Aj{ujr) = A\-{ur) + ^ifj(wr) (j = 1, ■ ■ ■ , 4). 

From the Equation (38), it is known that the amplitude of the velocity error can be expressed as 
a function of the inertial sensor error and the rotating angular rate. Therefore, for the purpose of 
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ensuring the minimum of the velocity error amplitude, the rotating angular rate should match the error 
of the inertial sensors. In order to verify and determine this matching relationship, we will design three 
groups of inertial sensor performances and calculate the optimal rotating angular rate of each group 
with MATLAB. Taking "MARINS", a high-precision INS manufactured by IXBlue, as an example, its 
positioning accuracy is 1 nm/24 h. It can be inferred that the constant bias and the scale factor error of its 
gyroscopes are lower than 0.01 °/h and 1 x 10~^ (°/h)^^, respectively. On the other hand, considering 
the real INS, the constant bias accelerometer is lower than 9.8 x 10^'' m/s^, and the installation error of 
gyroscopes is lower than 1 x 10"^ °/h. Therefore, on the basis of the above, the simulation parameters 
are set as Tables 1-3, respectively. 

Table 1. Inertial sensor performance (Group 1). 

Sensor Error Parameter Value 

Constant bias of gyroscope 0.01 (°/h) 

Constant bias of accelerometer 9.8 x 10~^ (m/s^) 

Scale factor error of gyroscope 1 x 10~^ (°/h)~^ 

Installation error of gyroscope 1 x 10~^ 



Table 2. Inertial sensor performance (Group 2). 

Sensor Error Parameter Value 

Constant bias of gyroscope 0.05 (°/h) 

Constant bias of accelerometer 4.9 x 10"^ (m/s^) 

Scale factor error of gyroscope 0.2 x 10"^ (°/h)~^ 
Installation error of gyroscope 5 x 10~^ (°/h) 



Table 3. Inertial sensor performance (Group 3). 

Sensor Error Parameter Value 

Constant bias of gyroscope 0.002 (°/h) 

Constant bias of accelerometer 1.96 x 10^"^ (m/s^) 

Scale factor error of gyroscope 5 x 10"^ (°/h)~^ 
Installation error of gyroscope 2 x 10^^° 



Figure 7 depicts the relationship between the amplitude of the velocity error and the rotating angular 
rate by MATLAB. The value of the horizontal axis corresponding to the minimum of the curve is the 
optimal rotating angular rate. The optimal rotating angular rates of these three groups are summarized 
in Table 4. 
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Figure 7. The relationship curves between the velocity errors and the rotating rates. 
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Table 4. Optimal rotating rates with different simulation conditions. 



Group 


Optimal rotating angular rate 


1 


21.6 7s 


2 


64.8 7s 


3 


7.2 7s 



Figure 8. IMU rotating scheme. 
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4. Numerical Simulations 

The optimal design method of the rotating angular rate is proposed in Section 3. The theoretical 
optimal rotating angular rates corresponding to three different groups of inertial sensor errors are 
obtained in terms of this method. In order to verify the correctness of the former conclusions, simulations 
are carried out in this section. The rotating scheme is proposed as shown in Figure 8. The turntable 
rotates the IMU back and forth along the azimuth axis through four orthogonal positions. The dwell 
time in each position is 45 s [25]. 

The performances of inertial sensors are described in Tables 1-3. The rotating angular rates are 
designed as 7.2 7s, 21.6 7s, 45 7s and 64.8 7s; as 21.6 7s, 64.8 7s and 7.2 7s are the optimal 
rotating rate of Groups 1-3, respectively, and 45 °/s is chosen for comparison. 

Figures 9-11 show the simulation results of the position error of different groups. The simulation 
duration is 72 h. 

Figure 9. Simulation results of Group 1. 
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Figure 10. Simulation results of Group 2. 
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Figure 11. Simulation results of Group 3. 
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Tables 5-7 summarize the simulation results. 



Table 5. Position error (Group 1). 



Rotating angular rate (°/s) 7.2 21.6 45 64.8 



Position error (nmile/72h) 1.2 0.6 1.5 1.3 
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Table 6. Position error (Group 2). 



Rotating angular rate (°/s) 7.2 21.6 45 64.8 
Position error (nmile/72h) 5.0 0.9 1.0 0.6 



Table 7. Position error (Group 3). 

Rotating angular rate (°/s) 7.2 21.6 45 64.8 
Position error (nmile/72h) 2.5 3.0 5.0 4.2 



It can be seen from Figure 9 and Table 5 that the position error of Group 1 with a rotating angular 
rate of 21.6 °/s is lower than the other rotating angular rates. From Figure 10 and Table 6, it can be 
seen that the position error of Group 2 with a rotating angular rate of 64.8 °/s is lower than the other 
rotating angular rates. From Figure 10 and Table 6, we can see that the position error of Group 3 with 
a rotating angular rate of 7.2 °/s is lower than the other rotating angular rates. The simulation results are 
in accordance with the conclusions of Section 3. 

5. Experimental Study 

In order to verify the effectiveness of the method, experiments were done in Songhuajiang River 
in Harbin, as shown in Figure 12. Experimental equipment, as shown in Figure 13, included SINS, 
single-axis turn-table, GPS, computers and a UPS. The SINS is manufactured by Harbin Engineering 
University, and the main parameters are shown in Tables 8 and 9. The single-axis turn-table is 
manufactured by Beijing Precision Engineering Institute for Aircraft Industry, and the main parameters 
are shown in Table 10. 



Figure 12. Experimental ship. 
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Figure 13. Single-axis RSINS. 




Table 8. The main parameters of the fiber optics gyroscope. 



Performance Parameters 

Dynamic range ±100 °/s 

Bias stability < 0.003 ° /h. 

Random walk < 0.001 VVh 
Scale factor error < 5 ppm 



Table 9. The main parameters of the accelerometer. 



Performance Parameters 

Dynamic range ±4 g 
Bias stability < 1 x 10"^ 



Table 10. Performances of the 920E single axis velocity and position turntable. 



Performance 


Parameters 


Diameter 


450 m 


Carrying capacity 


weight 50 kg 


Rotating accuracy 


±2" 


Rotating range 


continues and infinite 


Location accuracy 


±3" 


Location definition 


0.0001° 


Speed range 


0.005 - 200 °/s 


Speed accuracy 


5 X 10^5 (360°average) 


Speed accuracy 


5 X 10^"^ (10°average) 


Speed accuracy 


1 X 10~^ (l°average) 
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5.1. Experimental Scheme 

With the optimal design method proposed in this paper, the optimal angular rate for the experimental 
system is calculated to be 7 °/s. Therefore, three angular rates, 2 °/s, 7 °/s and 20 °/s, are used in the 
experiments for comparison. The experimental duration is 72 h, and the rotating scheme is the same as 
that in Section 4.1. 

The rotating angular rate is not constant in a real system, since the IMU needs to reverse the rotating 
direction. The positioning accuracy of RSINS is influenced by the IMU turned angle in the variable 
angular rate process [26]. Additionally, the analysis results show that the error will be minimum in the 
process of the varied angle rate when the IMU's turned angle achieves 0.27r. In order to eliminate the 
effects of the variable angular rate, the angular accelerations are set as 0.1 °/s^, 0.7 °/s^ and 5.6 °/s^, 
respectively, in the experiment. 

5.2. Experimental Results and Analysis 

Figures 14-16 show the position error of three experiments. 

Figure 14. Experimental results (I). 
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From the theoretical analysis of Section 3, it is known that the rotation modulation method can only 
reduce the amplitude of the oscillation error and the slope of the ramp error, but not eliminate them. 
Therefore, it can be seen from Figures 14-16 that the resulting position error comprises a ramp error 
with superimposed Schuler oscillation and Earth rotation oscillation, and the oscillation periods are 
84.4 min and 24 h. From Table 1 1, it can be seen obviously that the position error of Experiment II is 
much less than the other two. The experimental results are consistent with the theoretical analysis of 
Section 3. 
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Figure 16. Experimental results (III). 
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Table 11. Experimental results. 



Rotating angular rate (°/s) 2 7 20 



Position error (nmile/72h) 14.7 6.4 19.6 
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6. Conclusions and Future Work 

This paper presented a rotating angular rate optimal design method of the RSINS. Firstly, the RSENS 
and the principle of the rotation modulation were introduced, and the inertial sensor errors were modeled. 
Secondly, the relation between the velocity error and the angular rate was worked out using the Laplace 
transform and the inverse Laplace transform. The optimal rotating angular rate method was also proposed 
based on the relationship above. Thirdly, simulations were conducted, and the results confirmed that the 
optimal rotating angular rate could be computed correctly with this method. Navigation experiments 
were done, and the results of the positioning error confirmed the effectiveness and correctness of 
the theoretical analyze. However, the experimental time was just 72 h, because of the limitations of 
experimental conditions. As a result, we plan to do more further experiments, which will last a longer 
time, if the experimental conditions permit, in the future. On the other hand, the results of the simulations 
were superior to that of the experiments. The reason may be that an error of the initial alignment, which 
is ignored, in the simulations may occur. The initial alignment error is not the main factor of the rotating 
angular rate optimal design method from the analysis of simulations and experimental results. Even so, 
we will model the initial alignment error and focus on its impact on the RSINS in the future. 
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